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Abstract 

Critical properties of the Calogero-Sutherland model of BCN-type {BC jy-CS 
model) are studied. Using the asymptotic Bethe-ansatz spectrum of the BC^- 
CS model, we calculate finite-size corrections in the energy spectrum. Since 
the BCn-CS model does not possess translational invariance, the finite-size 
spectrum acquires the contributions coming from "boundaries" . We show that 
the low-energy critical behavior of the model is described by c = 1 boundary 
conformal field theory. Thus the universality class of the model is identified 
as a chiral Tomonaga-Luttinger liquid. 
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I. INTRODUCTION 



The Calogero-Sutherland (CS) models describe one- dimensional quantum many- 
body systems with inverse-square long-rang interactions. Among many variants of the CS 
model 0, a class of models which are not translationally invariant has been known over 
the passed years ||^. In particular the so-called CS model of BCj\f-tjpe (abbreviated as 
the BCn-CS model hereafter) is the most general model with interacting particles. The 
BCn-CS model is intimately related to the root system of type BCn and invariant under 
the action of the Weyl group of type Bj^j. Namely, the model is invariant under coordinate 
transformations 

(?i,92,---,giv) ^ (eig<x(i),e2g<x(2),---,eAf9<T(jv)), (1) 

where {qi,q2, ■ ■ ■ , ^at) G denote the coordinates of particles, ej E {±1} and a is an 
element of the symmetric group of A^ letters. Roughly speaking, the Weyl group of type Bjsi 
consists of the ordinary exchange of particle coordinates and the sign change of coordinates. 
As we will see below the latter is understood as the mirror image of particles with respect 
to a boundary. 

Recent works have made it clear that the BC ^-CS model is relevant to one-dimensional 
physics with boundaries. For instance, it was pointed out that the non-relativistic dynamics 
of quantum sine-Gordon solitons in the presence of a boundary is described by the BCn-CS 
model (with sinh-interaction) This model is interesting in view of the quantum electric 
transport in mesoscopic systems The Haldane-Shastry model, which is the discrete 

version of the CS model, with open boundary conditions can also be constructed by utilizing 
the root system of type BCn HjlOl- shall present further evidence for the relevance of 



the BCn-CS model to our understanding in one- dimensional physics including boundary 
effects. 

In this article we will analyze the long-distance critical properties of the BCn-CS model. 



Since the exact energy spectrum of the model is available [11|, we may apply the method of 
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finite-size scaling developed in conformal field theory (CFT) to study the critical behavior. 
The same technique has already been employed when the critical properties of the CS model 
of A^^i-type were considered ||12[. The universality class of the Aat-i-CS model is identified 
as a Tomonaga-Luttinger liquid which is equivalent to c = 1 Gaussian CFT. In what follows 
we will show that, in contrast to the Ajv-i-CS model, the BC^-CS model exhibits the critical 
behavior described by c = 1 CFT with boundaries [|l^ . Hence the universality class will be 
found to be a chiral Tomonaga-Luttinger liquid |jl4 . 



In the next section we first introduce the BC n-CS model and review the energy spectrum 
of the model obtained by using the asymptotic Bethe-ansatz. In section 3 we consider the 
thermodynamic properties. In section 4 the finite-size scaling analysis of the energy spectrum 
is performed. Finally, in section 5, we discuss various critical exponents of correlation 
functions. 

II. THE BCn-CS model 

Let us write down the Hamiltonian of the BCn-CS model 0. We put the system in finite 
geometry with linear size L and impose periodic boundary conditions. The Hamiltonian is 
then given by 

i=i"^i i<j<fc<7V [ sm -{qj-qk) sm —{qj + qk) 

/7r\2 ^ 1 /7r\2 ^ 1 

+X,{X, + 2X, - 1) - E -^jr- + 4A2(A2 - 1) y E ^1^' (2) 

j=ism'-qj i=isin^-2gj- 

where A, Ai and A2 are coupling constants which are assumed to be non-negative. It is 
clearly seen that the Hamiltonian is invariant under the action of the Weyl group 
of type B]\f. There exist several interaction terms which will need explanation. The term 
1/ sin^ {tt / L){qj + q^) expresses the two-body interaction between the j-th particle and the 
"mirror-image" (we place a mirror at the origin g = 0) of the k-th particle (j 7^ k). The term 
1/ sin^(7r/L)gy may be interpreted as the potential due to impurity located at the origin. 



The term 1/ siv? {ji / L)2qj describes the interaction between the j-th particle and its own 
"mirror-image" . All these terms required by invariance under the action of the Weyl group 
of type Bn violate translational invariance. Therefore, the total momentum is not a good 
quantum number for the BCn-G^ model. 

The Hamiltonian (0) can be cast into another form just by using the elementary identity 
sin 2 A = 2 sin A cos A. One gets 

^ 92 /7r\2 f 1 1 

w = -i:4 + 2A(A-i) Y. ; + -TT7 : 

j=i^'ij i<j<fc<7V [sm -{qj~qk) sm -{qj + qk) ^ 

^ /vr\2 JX, 1 ^ ^ /vr\2 JX, 1 

+ 7 + 7 E^^> (3) 

^-^"^ j=i sin —qj ^-^^ j=i cos —qj 

ij ij 

where /i = Ai + A2, v = A2. In this form of the Hamiltonian the term 1/ sin^^n/ L){qj + q^) 

is regarded as the boundary potential as before, while the last two terms in are regarded 

as the impurity potentials with the strength determined by and v respectively. The 

Hamiltonian is suitable for our present considerations. 

The eigenvalues of the Hamiltonian (j^) of the BC^-C^ model have been obtained by one 

of the authors The energy spectrum so obtained is shown to be reproduced exactly 



with the use of the asymptotic Bethe-ansatz (ABA) method [|TT|]. Let us recall the ABA 



formula for the -BCat-CS model. First of all the total energy of the system takes the form 

N 

en=y: (4) 

where pseudomomenta fcj's satisfy fci > A;2 > • ■ ■ > A^at > and obey the ABA equations 

N 

kjL = 2nlj + 7r(A — 1) ^ {sgn(fcj — ki) + sgn{kj + ki)} 

+7r{fi - l)sgn(/cj) + Tiiu - l)sgn(fcj), j = 1,---,N, (5) 



^Precisely speaking, this reference treated the case with u = (the Bj\[-CS model). However, we 



can easily obtain the formula for the BC n-CS model. The spectrum was also derived in [|lO|. 
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with sgn(x) = 1 for x > 0, = for x = and = —1 for x < 0. Here Ij {j = 1, ■ ■ ■ , N) 
are positive integers with Ji > J2 > ■ ■ ■ > /at > 0. These are quantum numbers which 
characterize the excited states. 

We emphasize here that, in contrast to the Atv-i-CS model, the Fermi surface of the 
BC N-CS model consists of a single point. This is due to the fact that pseudomomenta kj 
which are solutions to (H) are distributed only over the semi-infinite region as is shown in 
Fig. la. Therefore, in view of the bosonization picture, it implies that the low-energy critical 
behavior of the BCn-CS model will be effectively described by a left (or right)-moving sector 
of CFT (see Fig. lb). In addition to this, we also notice that the form of our Bethe-ansatz 
equations (^) is quite close to that appeared in the studies of the nonlinear Schrodinger 



equation on the half line |15,16| as well as the XXZ model with open boundary conditions 



|T^,[T8|. The critical behavior observed in these models p!5|-[18| is well described by boundary 



CFT [|T3|. It is inferred from these points that boundary CFT will play a role in our study 
of the BCj\f-CS model. 

Finally we rewrite our ABA equation (|^) for further convenience. As has already been 
mentioned, all the pseudomomenta kj are positive. However, one can make a trick so that kj 
takes values in (—00, 00) as in the bulk system. To realize this let us define J_j = —Ij, Iq = 
0, k-j = —kj and ko = with j = 1, ■ ■ ■ , A^, then we have 

1 1 ^ 

+^(/i + 1^- 2)sgn{k,) - |(A - l)sgn(2/c,) - ^(A - l)sgn(/c,), (6) 

where j = —N, —N + 1, ■ ■ ■ , N. The last two terms in (^ arise since the summation in (^ 
does not include the terms / = j and / = 0. Now the system turns out to have linear size 2L 
and the number of particles becomes 2A^ + 1. Note that the density of the system does not 
change. This doubling trick is known to be efficient when studying one- dimensional physics 



with boundaries [p!5|-[T8 
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III. THERMODYNAMIC PROPERTIES 



The purpose in this section is to discuss thermodynamics of the BC n-CS modeL Let us 
first consider the system at zero temperature. All the states inside of the interval [—kp, kp] 
are occupied, where the Fermi momentum kp is defined as kp = maxj/cj}. The thermody- 
namic limit is taken by 2L —>■ oo, 2N + 1 ^ oo with the density {2N + 1)/2L fixed. As 
usual we define the density of states by 

lim ^ , . =p{k), (7) 

Lh-»oo 2L[kj — Kj+ij 

and the sum is converted into integral 

— E - / dkp{k){ ). (8) 

ZIj J -kp 

From (H), (H) and ^sgn(x) = 25(x), it is shown that 

/kp 271 
dk'S{k - k')p{k') + — (/i + v- 2\)5{k), (9) 
-kp Li 

where the boundary effect manifests itself in the last term (~ ^/L). Notice that even for 
/i = 1/ = 0, it still modifies the equation. Upon taking the thermodynamic limit one can 
neglect the boundary term. The resulting equation is the same as for the Aat.i-CS model 
1^. Then it is immediate to get 

m = jij. (10) 

kp = 2Ti\d, (11) 
where we have put d = N/L. It is also straightforward to compute the ground-state energy, 

^ rkp 

= E i'^ff = 2^ / dkPpik) = 2L ■ e(o) (12) 

j=-N ■'-'^P 

with e(°) = An'^X^d^/'i in the 2L ^ oo limit. 

It is not difficult to extend the above analysis to the finite temperature case. At finite 
temperatures the pseudomomenta distribute over the infinite region (— oo, oo). One finds 
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ATi{p{k) + p^{k)) + 47r(A - 1) / dk'6{k - k')p{k' 



2n 

+ —{fi + u-2\)6{k), (13) 
where p^{k) is the hole density. Let 2L oo, then we have 

Following now the familiar procedure, we obtain the thermodynamic Bethe-ansatz equation, 

e{k) = e-p, + {X- l)Tlog |l + exp (-^e(A;)) | , (15) 

where T is the temperature, pc is the chemical potential and the energy density e{k) of 
particles is defined by 

'^^^ exp (-Uk)] . (16) 



p^{k) T 

Performing the low-temperature expansion of the free energy F(T) which is given by 

(F(T)-/i,(2iV+l))/(2L) = -— / dk\og{l + e-T<'')), (17) 

we have 

The second term in (^) is responsible for the linear specific heat C as T — 0. It is 
well recognized that the coefficient in C is universal modulo the Fermi velocity Vp which is 
not universal |]19|. In translationally invariant systems the Fermi velocity is determined by 



the dispersion relation. In the BC ^-CS model, however, one cannot rely on the dispersion 
relation since the momentum is not a good quantum number. So, in order to determine 
Vp, we have to take another point of view. As we observed, eqs.(|iy), (0) and coincide 



with those obtained in the Aat^i-CS model. Hence we may regard the y4jv_i-CS model as 
the bulk counterpart of the BCn-CS model. Since the Atv-i-CS model is described in terms 



of c = 1 CFT [n] we assume that the central charge for the BCn-CS model is also given 



by c = 1. Then, comparing C obtained from (|T8D to the formula C = 7TcT/{3vp) |T9| with 
c = 1 we find Vp = AnXd. We shall see in section 5 that the finite-size spectrum is in fact in 
accord with c = 1 CFT. 



IV. FINITE-SIZE SCALING ANALYSIS 



In this section we perform the finite-size scahng analysis of the energy spectrum of the 
BC N-CS model. To begin with, we summarize several fundamental formulas in boundary 



CFT 1 13 1 which we will need to analyze the energy spectrum. Let us first recapitulate the 
finite-size scaling form of the ground-state energy predicted by conformal invariance under 
free boundary conditions [|T^ 

= L6(°) + 2/ - , (19) 

where e*^"-* and / are, respectively, the bulk limits of the ground-state energy density and the 
boundary energy, v-p is the velocity of the elementary excitations. The Virasoro central charge 
c which specifies the universality class of the system appears as the universal amplitude of 
the 1/L term in (|19D. 

From the scaling behavior of the excitation energy one can read off the boundary critical 



exponents Xb [0]. This exponent Xb governs the power-law decay (parallel to the boundary 
surface) of a two-point function. Suppose a critical system on the half-plane {(?/,t) G 
R>o X R} with a surface at y = 0. {y is the perpendicular distance from a point {y, r) to the 
boundary and r means the imaginary time.) Let 0{y,T) be a local operator. We consider 
its two-point correlation function G{yi,y2,T) = {0{yi, Ti)0{y2, T2)) , which is a function of 
T = T1—T2 because of translational invariance along the surface. For |r| ^ yi, y2, we obtain 
the asymptotic form of G, 

G'(z/i,2/2,r)~4r- (20) 
To evaluate Xb we have to examine the scaling law 

E-E^'^ = ^Xb (21) 



with E being the excitation energy. It usually happens that the value of Xb is distinct from 
that of the bulk exponent for certain scaling operator. In terms of CFT, the bulk exponent 
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is expressed as the sum of left and right conformal weights, while the boundary exponent is 
equal to the left (or right) conformal weight. 

Let us now turn to the BC n-CS model. It is convenient to manipulate the ABA equations 
(^) directly. We can easily solve (^) to obtain 

2n 



[I^-(N~j + l)] + k 



(0) 

'i ' 



J = l 



(22) 



where 



,(0) 
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\{N-j) + 



(23) 



The ground state is thus specified by the quantum numbers if^ = N — j + 1, {j = 1, ■ ■ ■ , N), 
from which we get the Fermi point Ii'^ = N and the Fermi momentum kp = 2ti\N / L + 
7r(/i + u — 2\)/L. The ground-state energy is then obtained as 

N 



E 



(0) 
N 



E(*f) 

L 



^XN + 1a(/z + u - X)N' + ^ (3(/z + z/ - A)2 - A^) N' 



(24) 



.3 2 

We make a power expansion of ( pij ) with respect to 1/L while keeping the particle density 
d = N/ L fixed. The result reads 



E 



(0) 
N 



e^'^L + 2f + ^X{AN,f-^X, 



(25) 



where / = 7r^A(/i + u — X)d? and 



AN. 



fi + u — X 
2A ' 



(26) 



In (|25|) there appear no higher-order terms with L ™(m > 2). Note also the symmetric 
dependence of / and AN^ on /i, z/. 



There are several points which should be noticed in ([25|). First of all, besides the thermo- 
dynamic energy density e^^-* already computed in (|12|), one finds the boundary energy 2f in 
the term of order which is due to the absence of translational invariance in the system. 
The next order corrections proportional to 1/L turn out to provide valuable information on 
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"boundary effects". To see this, let us proceed a bit carefully by having decomposed the 
1/L-contributions into the last two terms in (pSf). We first recall that the size-dependence of 
the interaction is inevitably introduced for systems, as seen in (|^), when dealing with 
interacting particles in finite geometry. This gives rise to nonuniversal 1/L-corrections to 
the ground-state energy in addition to the universal one, as observed in the Aat.i-CS model 
1^ . In (|25|) , therefore, we think that the term —ttv^X/ (12L) suffers from such nonuniversal 



contaminations which, in direct comparison with ([19|) , yield the wrong value for the central 
charge. 

The other 1/L-correction term, 7rfpA(AA'^[,)^/L, is more interesting and understood as 
the "boundary effect" which consists of two kinds of contributions. As seen from (P), when 
we convert the BC^ system to the chiral system by using a trick of mirror image, we are left 
with particles moving only in one direction feeling the boundary potential depending on A, in 
addition to the impurity potential depending on /i and z/. These two types of scattering effects 
are combined into a quadratic form with respect to the "fractional quantum number" AN^ 
depending on both fi + u and A. Note that the quantum number ANf, physically represents 
the phase shift due to the scattering by the impurity- and boundary-potentials. Thus our 
ground-state energy E^^ is considered as the phase-shifted ground-state energy [^. If we 
imagine a hypothetical system which does not include these boundary contributions, the 
corresponding ground-state energy E^^^ is written as 

< = <'-^^(AiV.)^. (27) 

Having discussed the ground-state energy in detail, we next wish to calculate the finite- 
size corrections to the excited states. Looking at the ABA equations let us create an 
excited state by adding AA^ particles to the ground-state configuration. In this case, we 
have the pseudomomenta 

+ z/" 



\{N + AN-j) + 



(28) 



from which we immediately obtain the finite-size corrections to leading order in 1/L, 
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E 



(0) 



-C/AT — 



/i^AiV + (AiV + AiV,)2 - !I^A(AiV,)2, 



(29) 



where /x^*'-* = de'^^^ /dd = kp^ is the chemical potential. Note that this expression for the 
finite-size spectrum is essentially the same as that derived for the charge sector in the Kondo 
problem (see (49) in [0]). If we redefine e''^^ by e''^^ — /i^^^A^, we find 

27rt>p A 



E 



(0) 

N+AN 



E 



(0) 
N 



L 2 



(AA^ + ^N^Y . 



(30) 



Since any excitations which carry currents with large momentum transfer are barred due to 
the absence of translational invariance in the BC n-CS model, the remaining possible type 
of low-energy excitations are provided by particle-hole excitations labeled by non-negative 
integers n. The corresponding energy is simply obtained by adding 2TTVpn/L to (pO|). Hence 
we have 



E-E 



(0) 
N 



27rt>n 



^ {AN + ANbf + n 



(31) 



where E denotes the energy of the excited state specified by {AN, ANb,n). In the next 
section we argue that our result ( pT] ) is in accordance with the scaling law in c = 1 boundary 
CFT. 



V. CORRELATION FUNCTIONS 

Now that we have evaluated the finite-size corrections it is possible to read off various 
critical exponents using the scaling relation (^Tl). When comparing our result (^) with (^Tj) 



we have to replace L with 2L since L has been defined as the periodic length of the system. 
Bearing this in mind let us take an operator ipb which corresponds to the phase-shifted 



ground state. This operator can be assumed to be the boundary changing operator [20 



With this point of view, the phase-shifted ground state is an excited state relative to E^^^ 
in (p7|). The scaling dimension of ipb is obtained as 
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x,, = ^^{EP-EP) = l,iAN,)\ (32) 

where we have put ^ = ( = and hence AN^ = (^^ _ (^2y^2('^y 

We next consider an operator which induces the particle number change as well as the 
particle-hole excitation in the phase-shifted ground state. From ( ^OD and (|2l|) we have 

= ^ {e^nUn - eP) = ^ (AiV) ' + n, (33) 

where 

AiV = AA^ + ANb. (34) 

Scaling dimensions ( p2D and (|33|) take the form of conformal weights characteristic of c = 1 
CFT. The radius R of compactified c = 1 free boson is taken to be i? = ^. Let us concentrate 
on the self-dual point R = 1/ \pl {i.e. A = 2) where the symmetry is known to be enhanced to 
the level-1 SU{2) Kac-Moody algebra. In the BCn-CS model we have the other continuous 
parameters /i, u which should also be tuned to achieve the SU{2) point. It turns out that 
/i + 1/ = 0, 1, 2, 3 and 4 with A = 2 are the desired points. This follows from the following 
observations: When /i + z/ = 2 we have AA*";, = and hence 

X0 = ^(2AiV)2 + n (35) 

which is the conformal weight for the spin-0 irreducible representation of the level-1 SU (2) 
Kac-Moody algebra. When /i + z/ = 4 or we get ANi, = ±1/2 and thus 

= ^{2AN + if + n (36) 

which is the conformal weight of spin- 1/2 irreducible representation. When /i + z/ = 3 or 1 
we have ANb = ±1/4? thereby 

x^ = ^{AAN+lf + n. (37) 
This is the conformal weight for the unique irreducible representation of the level-1 twisted 



SU{2) Kac-Moody algebra The highest-weight state with x,^ = 1/16 is a twist field in 
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c = 1 CFT. Several SU{2) points identified in |T0| are in agreement with our result. Thus we 



conclude that the low-energy critical behavior of the BC n-CS model is described in terms 
of c = 1 boundary CFT, i.e. the universality class of a chiral Tomonaga-Luttinger liquid. 

Further considerations on the low-energy critical properties of the BC n-CS model require 
a clear distinction between two pictures corresponding to two possible sets of quantum 
numbers. One is a set of quantum numbers (AA^, AA*";,, n) and the other is a set of (AA^, n) 
where AA^ is regarded as the ordinary particle number change in (^) (forgetting about 
ANh in (|3^) ). The picture based on the set {AN, ANb,n) is relevant when describing the 
long-time asymptotic behavior of the system in which we suddenly turn on the boundary 
effects in the ground state. The X-ray absorption singularity in the Kondo problem, for 
instance, is considered in this type of picture p0|j21|| . The boundary changing operator iph 
is described in this picture with {AN, ANi,,n) = {0, ANi,,0). If we use the set (AA^, n) 
instead, our picture is independent of ( and adequate to compute the critical exponents of 
ordinary correlation functions with boundary effects. 

Let us consider the one-particle Green function in the above two pictures. Let 
(AA^, ANh, n) = (1, AA";,, 0) in the first picture. This choice of quantum numbers determines 
the long-time asymptotic behavior of the field correlator (the one-particle Green function) 
when boundary potentials are turned on at r = 0, 

(vI/t(r)M/(0))sMdcn ~ (38) 

where 

^G = ^{l + Am)' = ^(l + Q . (39) 



Here (■ ■ ■)suddcn stands for the expectation value when the boundary potential is suddenly 
switched on. On the other hand, if we let {AN,n) = (1, 0) in the second picture, the field 
correlator takes the form, 

(vl/t(r)vl/(0)) ~ (40) 
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where 



^3 = (41) 



which describes the ordinary one-particle Green function. In this case, the boundary critical 
exponent Xg linearly depends on A. Contrary to these Green functions, the density- density 
correlation function is controlled by the excitations which do not change the number of 
particles. Hence, it should have the long-time asymptotic form, 

(p(t)p(O)) ~ (42) 

which follows by taking the quantum number {AN,n) = (0, 1) in (33). Note that there do 
not appear anomalous exponents in this correlator. One can easily see that this is also the 
case for sub-leading terms r^^'^ in which the quantum number is chosen as (AA^, n) = (0, k). 
This fact will be confirmed shortly in the following. 

We now compare our result with the explicit calculations of the dynamical correlation 
function. In the case A = 1, z/ = with n arbitrary which corresponds to the noninteracting 
system, the dynamical density-density correlation function for the BCn-CS model has been 
obtained by Macedo ||23| (see also [^]). In the thermodynamic limit, the density-density 



correlation function G{yi,y2,T) has the form 

TT^ f°° _i 2 

G{yi,y2,T) = —ym duie 2^ (tti/iv/mT) (vryav^) 

/■I 1 2 

x/ rfM2e2^™V i(7r?/iv/M^)J i(7r?/2Vw2), (43) 

where Jv{z) is the Bessel function and r is the imaginary time. When fi = 1/2 + m [m = 
0, 1, ■ ■ ■) it is not difficult to evaluate the large-r asymptotic behavior by making use of the 
series expansion of Jm{z)- After some algebra we obtain 

CO /i\2fc 00 /1\ '+"^+2 

G(2/i,2/2,r) = EA(z/.) - +T.Bi{y,)[-] e-^-'^ (44) 



k=l ^'^ ^ 1=0 ^'^ 



where Ak{yj), Bi{yj) are some functions. As r — > 00 with yi, y2 fixed, the second term 
vanishes exponentially, yielding 
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+ $ + $ + (45) 

Notice that the exponents are independent of m {i.e., /i). The density-density correlation 
function is considered in the picture based on (AA^, n). Then we see from (|33|) that all 
these exponents are precisely understood in terms of the excitations (AA^, n) = (0, k). This 
means that the correlation function G is dominated by the particle-hole excitations, and 
hence there is no way of depending on A. Therefore the result (^) completely agrees with 
our prediction by CFT analysis. We are thus led to conclude that the power-law decay in 
(^) is universal irrespective of A (but with z/ = fixed) though ( ^Sj) is verified at A = 1. We 
stress that this remarkable feature in the density-density correlation function is inherent in 
chiral Tomonaga-Luttinger liquids |1l4| . 

Finally we briefly mention possible applications to the (chiral) random matrix theory 
| 25| . Let us recall the Bn Calogero-Moser model {Bn-CM model) in the rational form ||^, 



^ r 1 1 

^c-M = -E7^ + 2A(A-l) E 



j=i i<j<k<N I ^''^)^ ("^i •'^'^y 



N N 

-f/i(/i-l)E- + ^'E^?' (46) 
j=l j=l 

with uj > 0. In the thermodynamic limit, this model belongs to the same universality class 
as the Bjy-CS model which is equivalent to the BCjy-CS model at u = 0. The ground-state 
wave function for the B^-CM model takes the form of Jastrow-type 

^ / 1 \ 

¥'\x,,X2,---,xn)=^ n k?-4l^ni^'l^exp (47) 

i<j<k<N 1=1 \ I J 

where A/" is a calculable normalization constant. Notice that \E'(°)(a;i, X2, ■ ■ ■ , a;Ar) de- 
pends only on the a;|'s. Then, introducing new variables Zj = x^, one should note that 
|^'^°^(xi, X2, ■ ■ ■ , xn)]"^ is identical to the probability distribution function for the eigenvalues 
Zj of the Laguerre ensemble when A = 1/2, 1 and 2 (with appropriate values of fi and u) cor- 
responding to the ensembles of orthogonal, unitary and symplectic types [^, respectively. 
Therefore, it will be very interesting if the long-time asymptotic behavior of correlation 
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functions in the Bn-CM model obtained in the present work is directly compared with the 
results in the Laguerre random matrix theory. 

In summary, we have investigated boundary critical phenomena in the BCn-CS model. 
The boundary effects come from both the impurity potentials and interactions between 
particles and "image" particles. Making use of boundary CFT, we have obtained boundary 
critical exponents, and clarified the critical properties of the BCn-CS model in terms of 
chiral Tomonaga-Luttinger liquids. 
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FIGURES 





FIG. 1. (a) The Fermi surface consists of a single point k = kp. (b) Schematic illustration for 
the bosonization picture. 
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